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We present a model for theoretical description of the dynamics of a system of spinor cavity 
polaritons in real space and time accounting for all relevant types of the interactions and effective 
magnetic fields. We apply our general formalism for the consideration of the polarization dynamics of 
the coherently driven one dimensional polariton channel. We investigate effect of the temperature 
and the longitudinal-transverse splitting on the spin (polarization) multistability and hysteresis 
arising from the polarization-dependent polariton-polariton interaction. We show that the effect 
of the phase of the driving laser pump is as important as its strength and demonstrate that the 
multistability behavior can survive up to high temperatures in presence of longitudinal-transverse 
splitting. 



I. INTRODUCTION 



Cavity polaritons are composite particles, arising from 
strong coupling between photonic mode of a planar semi- 
conductor microcavity and exciton transition in Quan- 
tum Well (QW), embedded in the cavity at the point 
where the electric field of the confined electromagnetic 
state reaches its maximum. Having a hybrid half-light - 
half-matter nature cavity polaritons demonstrate a set of 
peculiar properties which make them different from other 
quasiparticles in solid state systems. Extremely low ef- 
fective mass of the cavity polaritons (about 10~ 4 — 10 -5 
of the free electron mass) together with strong polariton- 
polariton interactions makes polaritonic system an ideal 
candidate for observation of the variety of quantum col- 
lective phenomena at surprisingly high temperatures. 
Achievement of polariton Bose-Einstein condensation 
was first reported at T = 20A' 1 and later on even at room 
temperature^. Later on polariton superfluidity 3 ., Joseph- 
son effect and formation of topological excitations 5 
were experimentally observed. Other theoretically pre- 
dicted effects such as polariton self-trapping^, polariton- 
mediated superconductivity^ still wait for their experi- 
mental confirmation. 

Besides fundamental interest, quantum microcavitics 
in the strong coupling regime can be used for optoelec- 
tronic applications^. For more then a decade, the only 
object of the study in this context was polariton laser^- 
a novel type of the coherent emitter which explores the 
possibility of the polariton BEC. In the last years, how- 
ever, the emphasis became to shift on other types of the 
devices based on transport properties of cavity polari- 
tons in real space. It was noticed that the peculiar spin 
structure of polaritons opens a way for creation of optical 
analogs of spintronic components (so-called spinoptronic 
devices^). With respect to optics, spin-optronics has the 
advantage of being able to use particle-particle interac- 
tions occurring in nanostructures and resulting in strong 
nonlinearities. With respect to spintronics, it has the 
advantage of strongly reducing the dramatic impact of 
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FIG. 1: Sketch of the system showing the position of the 
polariton channel with respect to the Poincare sphere (also 
known as Bloch sphere and pseudospin sphere). The latter 
serves to illustrate different possible polarizations of the light. 
If the vector of the pseudospin lies in the xy plane, the light is 
linearly polarized and if it is parallel to z-axes it is circularly 
polarized. Other orientations correspond to the general case 
of elliptical polarization. Red arrow shows the direction of the 
effective magnetic field created by the longitudinal-transverse 
splitting f2. 



carrier spin relaxation or decoherence, which has severely 
limited the achievement or the functionality of any work- 
ing semiconductor-based spintronic devices. 

In this context, the analysis of one-dimensional (ID) 
polariton transport is of particular importance^, as ID 
polariton channels are fundamental building blocks of 
such future spinoptronic devices as polariton neurons^ 
and polariton integrated circuits^ 3 -. It should be noted, 
that current state of growth technology offers a large 
variety of methods of the lateral confinement of cavity 
polaritons^ and polariton quantum wires (ID polariton 
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channels) can be routinely produced. 

Currently, the theoretical study of transport of spinor 
cavity polaritons in the real space is based on the as- 
sumption of full coherence of the polaritonic system. 
Polariton-polariton interactions are either neglected^—, 
either treated within frameworks of spinor Gross- 
Pitaevskii equations (GPe)^. The non-coherent pro- 
cesses coming from the interaction of the polaritonic sys- 
tem with a phonon bath in most cases are not accounted 
for or treated within simple phenomenological models 
lacking microscopic justification 1 ^— . On the other hand, 
it is clear that polariton-phonon interaction is of crucial 
importance, as it provides a thermalization mechanism 
for a polaritonic ensemble and drastically affects such 
experimentally observable quantities as first and second 
order coherencie s 22 ' 23 . 

It should be noted, that for spatially homogeneous po- 
lariton system polariton-phonon interactions can be ac- 
counted for using a system of the semiclassical Boltz- 
mann equations 2 ^ - — . This method, however, has several 
serious drawbacks. First, it is based on the assumption 
that the system is fully incoherent, and thus variety of 
intriguing nonlinear phenomena such as bistability and 
multistability can not be described. Second, it provides 
the information about occupation numbers in the recip- 
rocal space only, and thus can not be used for description 
of the dynamics of the spatially inhomogeneous system. 
This makes this formalism inappropriate for modeling of 
spinoptronic devices based on polariton transport in real 
space. 

In the present paper we present formalism suitable for 
the description of the dynamics of an inhomogeneous 
spinor polariton system in real space and time accounting 
for all relevant types of the processes. Namely, we take 
into account polariton-polariton and polariton-phonon 
interactions and effective longitudinal-transverse (TE- 
TM) magnetic field acting on polariton spin. Our con- 
sideration is based on the Lindblad approach for density 
matrix dynamics and represents a generalization of our 
previous work where spinless case was considere d 28 ' 29 . 
We use our results for modeling of the spin dynamics 
of the polaritons in ID channels, investigating the role 
played by dccohercncc at different temperatures. 



II. FORMALISM 

We describe the state of the system (polaritons plus 
phonons) by its density matrix for which we apply 
the Born approximation factorizing it into the phonon 
part which is supposed to be time-independent and cor- 
responds to the thermal distribution of acoustic phonons 

X P h = exp | — | and the polariton part Xpol whose 

time dependence should be determined, \ = Xph ® Xpol ■ 
Our aim is to find dynamic equations for the time evo- 
lution of the single-particle polariton density matrix in 



real space and time 

p^(v,v',t) = Tr{^t( r>t )^,( r ' jt ) p J (l) 

= {ii{v,t)^(v',t)) 

where ^(r, t), Vv ( r > t) are operators of the spinor po- 
lariton field, the subscripts a, a 1 = ±1 denote the z- 
projection of the spin of cavity polaritons and correspond 
to right-and left-circular polarized states and the trace 
is performed by all the degrees of freedom of the sys- 
tem. The particularly interesting quantities are matrix 
elements with r = r' which give the density and polar- 
ization of the polariton field in real space and time 

n(r,t) = ^ P<y,a{r,r,t), (2) 

cr=±l 

s z{r,t) = i [p +h+1 (r,r,t) - p-x _i(r,r,i)] , (3) 
s x (r,t) + is y (r,t) = p +1 _i(r,r,i). (4) 

The off-diagonal matrix elements with r/r' also have 
physical meaning and describe spatial coherence in the 
system. 

To obtain expressions for temporal dynamics of the 
components of single particle density matrix it is con- 
venient to go to the reciprocal space, making a Fourier 
transform of the one-particle density matrix, 

p CT , CT ,(k,k',t) = (2ir) d /L d J e^ kr - k ' r ')p CT)CT ,(r J r',t)drdr' 
= Tr { a ^k a ^.k'x} = (a£ k ( V,k'), (5) 

where d is the dimensionality of the system (d = 2 for 
non-confined polaritons, d = 1 for the polariton channel), 
L is its linear size, a^" k , a CT k are creation and annihilation 
operators of the polaritons with circular polarization a 
and momentum k. Note, that we have chosen the pref- 
actor in a Fourier transform in such a way, that the val- 
ues of p(k, k',t) arc dimcnsionless, and diagonal matrix 
elements give occupation numbers of the states in dis- 
cretized reciprocal space. Knowing the density matrix 
in reciprocal space, we can find the density matrix in 
real space straightforwardly applying the inverse Fourier 
transform. 

The total Hamiltonian of the system can be repre- 
sented as a sum of two parts, 

H = n 1 + n 2 , (6) 

where the first term T-L\ describes the "coherent" part 
of the evolution, corresponding to free polariton propa- 
gation, polariton-polariton interactions and the effect of 
TE-TM splitting, and the second term H2 corresponds 
to the dissipative interaction with acoustic phonons. The 
two terms affect the polariton density matrix in a quali- 
tatively different way. 
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A. Polariton-polariton interactions 

The part of the evolution corresponding to "Hi is given 
by the following expression 

kcr k,(T 

+ U l 2J °ki,cr k 3 ,cr a ki+P,ff k a -p,a 
ki,k2,p,(T 

+ C/2 2J a ki,o- a k 2 ,-o- a ki+P,o-ak 2 -p,-cr, 
ki,k2,p,(T 

where _Ek gives the dispersion of the polaritons, f2(k) 
is the TE-TM splitting corresponding to the in- plane 
effective magnetic field leading to the rotation of the 
pscudospin of cavity polaritons, Ui is the matrix ele- 
ment of the interaction between polaritons of the same 



circular polarization, t/2-matrix element of the interac- 
tion between polaritons of opposite circular polarizations. 
In the current paper we neglect the p-dcpendcncc of 
the polariton-polariton interaction constant coming from 
Hopficld coefficients for simplicity. As well, we will sup- 
pose Q (k) = const which corresponds well to the situa- 
tion of the polariton channel (but not for 2D polariton 
system). 

The effect oiTLi on the evolution of the density matrix 
is described by the Liouvillc-von Neumann equation, 

ih(d tX ) {1) = [Hr,x], (8) 

which after the use of the mean field approximation leads 
to the following dynamic equations for the elements of the 
single-particle density matrix in the reciprocal space (the 
derivation is completely analogical to those presented in 
Ref.il): 



-ih{d tPa , a , (k, k')} (1) = (E u - £ k ' W< (k, k') + f2 [p_ CT)CT , (k, k') - p a ,_ a , (k, k')] (9) 
+ Ui ^ [/w(ki,ki - p)p CTjCT /(k - p,k',) - p a ,^> (ki , ki - p)p CT>0 -< (k, k' + p)] 

ki,p 

+ U2^ [P-<T,-<x( k l7 k l - P)/W'( k - P> k ') - P-<r',- CT '( k l> k l - P)P<T, CT '( k 7 k ' + P)] ' 

ki,p 
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B. Scattering with acoustic phonons 

Polariton-phonon scattering corresponds to the inter- 
action of the quantum polariton system with the classical 
phonon reservoir. It is of dissipative nature, and thus 
straightforward application of the Liouville-von Neu- 
mann equation is impossible. One should rather use 
the approach based on the Lindblad formalism, which 
is standard in quantum optics and results in the master 
equation for the full density matrix of the system^. For 
the convenience of the reader, we give the main steps of 
the derivation of the dissipative part of dynamic equa- 
tions for spinor polariton system, omitting however all 
technical details which can be found elsewhere^. 

The Hamiltonian of the interaction of polaritons with 
acoustic phonons in Dirac picture can be represented as 

H 2 {t) = n-(t) + H + (t)= (10) 
= ]T ^(q)e l(Bk+ "- Bk)t a+ k+q ^,k(6 qe -^ t + b±^*), 

<T,k,q 

where a^k are operators for spinor polaritons, 6 q opera- 
tors for spinless phonons, E^ and w q are dispersion re- 
lations for polaritons and acoustic phonons respectively, 
D(q) is the polariton-phonon coupling constant. In the 
last equality we separated the terms H + where a phonon 
is created, containing the operators b + , from the terms 



7i~ in which it is destroyed, containing operators b. 

Now, one can consider a hypothetical situation when 
polariton-polariton interactions are absent, and the re- 
distribution of the polaritons in reciprocal space is due 
to the scattering with acoustic phonons only. One 
can rewrite the Liouville-von Neumann equation in an 
integro-differential form and apply the so called Marko- 
vian approximation, corresponding to the situation of 
fast phase memory loss (see Ref. |30| for the details and 
discussion of limits of validity of the approximation) 

(dtxf ] = ~ f dt> l<H 2 (t); [H 2 (ty, X (t)]) (11) 

= s AE [2 {n+xH- + n- x n+) 

- (n + n- + h-% + ) x-x (u + u- + %-%+)} , 

where the coefficient Sae denotes energy conservation 
and has dimensionality of inverse energy and in the cal- 
culation taken to be equal to the broadening of the po- 
lariton stated. For time evolution of the mean value of 
any arbitrary operator (A) = Tr(xA) due to scattering 
with phonons one thus has (derivation of this formula is 
represented in Ref. [H) : 

{d t (A)} (2) = 5 AE (([H-;{A;H+}}) + ([H+;IA- : H-}})) . 

(12) 
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Putting A = \a a i y in this equation we get the con- 
tributions to the dynamic equations for the elements of 



the single-particle density matrix coming from polariton- 
phonon interaction: 



(2) 



E 2VK(q) [(n k , CT + l)n k+qi(7 (n^ + 1) - n k , CT (n k+q>CT + l)ng 

- i(/3o. ,- ff (k, k)p_ CT)CT (k + q, k + q) + /9- CT . CT (k, k)p CT ,_ CT (k + q, k + q)) 

^ 2V^(q) [(n k , CT + l)n k+qj(7 n?.q - n k ^(n k+q , a + l)( ? i p ' q + 1) 

- \{p a ,- ff (k, k)p_ CT)CT (k + q, k + q) + /0- CT . CT (k, k)p CT _ ff (k + q, k + q)) 



(13) 



{9 tPff ,_ CT (k,k)}( 2 ) 

q!^k<Sk+q 

E 2W ^) 

q,-Ek>Sk+q 



/^-o-fk, k) ( -(n k+q!(T +nk+ q ,o-') 



p CT ,-er(k, k) ( -(n k+q>CT + n k +q,<r') , 



1 1' 1 ' 
'q 



— cr (k + q, k + q) ( -(n kjCr + n. k)(T /) - n 



(14) 







ph 


) 







{S t/ V, CT (k, k')} (2) = E K-( k ' k ')K +q .a - <) + P- CT , CT (k, k> CT ._ a (k + q, k + q)] 

- E ^(q) [/3 CT .a(k,k')(n k+q!CT + ri q ' l + l)+ / 9_ CT . CT (k,k')pa,- CT (k + q,k + q)] 

+ E W(q) [p v (k,k')(%^ q , ^ <) + p.-^k.k'lp^^k' + q,k' I q)] 

E W(q) [ P aA^ k')(n k ' +q , CT + nf + 1) + p CT ,_ CT (k, k')/9- CT , CT (k' + q, k' + q)] , 
q>s k />-E k ' + q 



(15) 



{ftftr.-^k, k')} (2) = E ^(l) [P--( k ' k ')K + q,. - O + P-a,-«r(k, ; k>„,_a (k + q, k + q)] (16) 

q^k<^k+q 

- E W^(q) [p< T ,_ CT (k,k')K+ q , ( x + n q ' l + l)+p_ (j! _ (T (k,k')/3 CT .- (J (k + q,k + q)] 

q^k>^k+q 

+ E [^.— ( k > k ')K'+ q , CT ' - <) + PaAK k> CT) _ CT (k' + q, k' + q)] 

q.B k /<B k / +q 

E W ^ [p°,-°(^ k')(n k / +q , CT ' + nf + 1) + p„ i<r (k, k> CT ,_ CT (k' + q, k' + q)] , 

q,B k ,>E k , +q 



r 



where, n k . CT = ^.^(k, k), j o + i : _i(k) = s x (k) + isy(k) in the reciprocal space which were obtained earlier us- 
and W(q) denote spin-independent scattering rates with ing another techniques^ 2 -. The equations for off-diagonal 



acoustic phonons (see Ref. |28| for the details). The first matrix elements with k ^ k' describe their decay which 

two equations corresponding to k = k' are nothing but physically corresponds to the decay of the coherence in 

the spinor Boltzmann equations for polariton-phonon the system coming from polariton-phonon interactions, 
scattering describing the redistribution of the polaritons 
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C. Pumping terms 

In this paper we concentrate on the case when system 
is pumped by external coherent laser beam. The corre- 
sponding Hamiltonian can be introduced as 

n cp = Y,p^A t )^ + h - c - ( 17 ) 

k.cr 

Here p k is the Fourier transform of the pumping ampli- 
tude in real space 

Mx,f) = P« r (x)e ik » x e-*-» t , (18) 
I 



D. Dynamics of the polarization 

The dynamics of the circular polarization degree p c of 
the light emission from the ground k = state can be 
defined as 



n k=0 n k=0 

where n k=0 (t) and n^ =0 (t) stand for the populations of 
polaritons with pseudospin ±1 correspondingly in the 
ground state of the dispersion. 

One should mention the effect of longitudinal- 
transverse splitting f2 on the polarization degree dynam- 
ics, since it couples a + and a~ modes together. Its role 
becomes more evident if one switches to the pseudospin 
formalism which is described in the introductory part of 
the manuscript. From the formal point of view, TE-TM 
splitting is equivalent to the effective magnetic field in 
+x direction Jl = e x fl (along the quantum wire). In the 
same time, the polariton-polariton interaction gives rise 
to another effective magnetic field oriented in z direction 
(structu re g rowth axes) fi p _ p cx e z {U\ — U2){n + — n~) 
(see Ref . I33T) . Therefore, the total effective magnetic field 
represents superposition l~2 fot = f2 p _ p + ft. Accordingly, 
it is possible to rewrite the kinetic equations as a cou- 
pled equations for occupation number n a and in-plane 
pseudospin Sj_— . Considering only the effect of effective 
magnetic fields (assuming infinite lifetime and absence of 



where -P<r(x) is the pumping spot profile in real space, k p 
is an in-plane pumping vector resulting from the inclina- 
tion of the laser beam as respect to the vertical, and w p 
is the pumping frequency of the single-mode laser. Time 
evolution of the arbitrary element of density matrix is 
given by2£ 



{cW(k,k')} (cp) - \{pl a {t){*v^)-PK^{t){a Ka y), 

(19) 

where the time evolution of the mean values of the anni- 
hilation operator reads 



(20) 



interaction with phonons), coupled equations are given 
as 

d t n+ =Q cx e z ■ {S ± x «), (22) 
d t S± cx (S± x ft p _ P ) + ^«= - « k =o)°- ( 23 ) 

This corresponds to the precession of the pseudospin 
along time dependent magnetic field, which leads to its 
non-trivial dynamics. 



III. RESULTS AND DISCUSSION 

We consider a microcavity based on AlGaAs family of 
alloys and use the following parameters. The Rabi split- 
ting was taken equal to 15 meV, polariton effective mass 
3xl0 -4 of the free electron mass and detuning between 
the pure photonic and excitonic modes 3 meV. The po- 
laritonic quantum wire is 50 fiia long and 2 fim wide. 
Further, we use typical polariton lifetime in a medium 
Q-factor microcavity, r = 2 ps. The polariton-polariton 
and polariton-phonon scattering rates have been taken 
independent on the wavevector for simplicity. The matrix 
element of polariton-polariton interaction was estimated 
using expression U ~ 3Ei,a^/S, where Eb is the exciton 
binding energy, as is its Bohr radius and S is the area of 
the wire, which gives U ~ 20 neV. The polariton-phonon 
scattering rate W = 10 s s _1 . Pump laser is detuned 



dt(ak.a) = —-zPk,o{t) - -E k (a k a ) - -Ok(ak,-o-) 
h n n 

- t. X](^i/V<x( k 2,k 2 - p) + C/2/0-<T,- o -(k 2 ,k2 - p))(a k + p , CT ) 

+ Yj W ( c [(«k+q,<r-^ ,l )(a k>CT )+p_ (rjCT (k + q,k + q)(a k _ <T )] 

q,-E k <i5k+q 

- VF(q)[(n k+q ^+riP /l + l)(a k ^)+ / 9_ CT . (T (k + q,k + q)(a k ^ (T )]. 

q,^k>^k+q 

I 
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FIG. 2: Dependence of the circular polarization degree of the 
driven polariton mode on the circular polarization degree of 
the driving pump in the absence of TE-TM splitting for dif- 
ferent temperatures. Hysteresis loops, the signatures of the 
multistability behavior, shrink with the increase of tempera- 
ture (T=0, 20, 60 K) and finally disappear at T^lOO K. Inset: 
dependence of polariton population versus pumping intensity 
for single component system demonstrating the phenomenon 
of bistability (T=0 K). 



above the energy of the lower polariton branch by 5 = 1 
mcV. We consider the case of spatially homogeneous cw 
pump of different polarizations. 

The bistable (for spinless condensate) and multistablc 
(if one accounts for the spin) behavior of a polariton sys- 
tem in ID and 2D quantum systems has already been 
investigated theoretically in a number of works (see, for 
instance, Ref. |34| ) and was reported experimentally^— . 
Most of the theoretical approaches are based on solution 
of the GPe. Unfortunately, this technique does not allow 
to account for the dissipation dynamics of polaritons due 
to interaction with the crystal lattice (phonon-mcdiatcd 
processes). The density- matrix approach, which is being 
developed in current manuscript, does. In the limiting 
case of zero temperature we immediately reproduce the 
results obtained by the GPe, as expected. 

The multistability (with multi-hysteresis) characteris- 
tic is shown in Fig. 2 for different temperatures in range 
1-100 K, in the absence of TE-TM splitting. In this case, 
since there is no mechanism of the transition between tr + 
and a~ modes, this effect can be understood in terms 
of independent bistable dynamics of a + and cr~ modes. 
Accordingly, in inset we present the population hystere- 
sis curve of a spinless polaritons to clarify the forthcom- 
ing discussion. Let's begin with inset. In the certain 
range of pumps the polariton population can take two 
different values depending on the history of the pumping 
process. If we slowly increase the intensity of pump, at 



some threshold value P^ h the population of the ground 
state jumps up abruptly due to the resonance of the 
blue shifted polariton energy with the energy of the laser 
mode. The system keeps staying at this high-populated 
state with further increase of the pump intensity. In the 
backward direction, when we decrease the intensity of 
pump, the bistable transition to the low-populated state 
appears at the lower pump intensity (p[^ < Pth) an d 
hysteresis curve appears. 

With account for spin, polariton-polariton interaction 
becomes polarization-dependent, which leads to multi- 
stability of the polariton circular polarization (see Ref. 
[HI for the detailed discussion of the situation at T = 0). 
This effect is illustrated in main plot of Fig. 2 for different 
temperatures, where the pump intensity is fixed and its 
circular polarization degree p p is being changed. Let us 
explain this phenomenon with the help of above discus- 
sion for spinless case. Keeping the total pump intensity, 
lets change its circular polarization from a~ (<f> = tt, and 
p p = —1) to a + (cf> — 0, and p p = 1), see Fig.QJ. Initially 
there exists only a~ polaritons in the ground state, thus 
p c = — 1. As p p is increased, a + component starts to 
become more populated, and at certain threshold value 
of p p the first bistable jump up in p c occurs that im- 
plies the abrupt increase of the <j + component. Further 
increase of p p leads to the second jump up of the polar- 
ization degree p c due to the bistability jump down of a~ 
component from high population state to low population 
state. Finally the system reaches the state with only a + 
component and p c — 1. In the backward direction (de- 
crease pump polarization degree from +1 to -1) the first 
jump down is due to the abrupt increase of o~ compo- 
nent, while the second jump down is explained by the 
abrupt decrease of a + component occupancy 

With increasing temperature, the multistability loops 
start to shrink and become totally destroyed at about 
T w 100 K. It occurs due to the dissipation pro- 
cesses coming from interaction with acoustic phonons. 
At higher temperature the spin-independent polariton- 
phonon interaction makes the dependence p c {pp) quasi- 
linear, as it should be expected, indeed, in the case when 
coherent nonlincaritics play no role and there is no tran- 
sition between circular polarized components the polar- 
ization degree of the system should coincide with those 
of the pump. 

Now let us introduce the TE-TM splitting to see its 
effect on polarization multistability. The corresponding 
term removes the isotropy in the xy plane since it acts as 
an effective magnetic field in +x direction. Consequently, 
the population of each component ri^ becomes dependent 
not only on circular polarization of the pump, but also 
on its in-plane component as can be seen from Eq. (|22[) . 
In Fig. 3, the dependence of the internal circular polar- 
ization degree of the system p c is plotted as a function of 
the circular polarization degree of the pump p p for three 
different in-plane angles 8 between the in- plane pseu- 
dospin of the pump and direction along the wire: axis 
Ox (see Fig. [T]). Azimuthal angle 8 comes in the pump- 
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FIG. 3: Internal versus external circular polarization degree 
for different xy-plane projections of the pseudospin of the 
pump (different azimuthal angles 8) for Q, = 0.08 meV, and 
T=0 K. Due to the finite value of the TE-TM coupling, the 
equivalency between x and y linear polarizations is broken. 
Therefore, the choice of meridian line of Poincare sphere along 
which the pump laser evolves from <r~ to a + state becomes 
crucial for the polarization dynamics of the polariton system. 



ing Hamiltonian as the relative phase factor between the 
pumping amplitudes, i. e. p+ = e l9 p-. It is observed 
that relative phase drastically modifies the profile of the 
Pc(pp) plot. First, we note that a finite y-component of 
the pump pseudospin (0 ^ 0, it) destroys the symmetry 
of the multistability curve with respect to the change of 
the sign of circular polarization of the pump (see dot- 
ted green line). Also, comparing the results for 8 = 
(+x direction) and = it (-x direction) cases, we see 
two quite different polarization behaviors. This differ- 
ence can be understood from the first term of the kinetic 
equation (|23|) for Sj_, where the two cases (8 = and 
8 = it) give contributions with opposite signs. Therefore, 
the internal circular polarization degree becomes highly 
sensitive to the choice of meridian on the surface of the 
Poincare sphere along which the circular polarization of 
the laser is changed between a + and o~ . It should be 
noted, that for circularly polarized pump the effect of ro- 
tation of pseudospin due to TE-TM splitting is invisible: 
for the used value f2 = 0.08 meV and circular polarized 
pump the effect of macroscopic self-trapping plays major 
role^^r— and rotation of pseudospin is blocked. For re- 
duced values of circular polarization self-trapping regime 
is lost and effects of TE-TM splitting become visible. 

Let us finally analyze the combined effect of the TE- 
TM splitting and scattering on phonons (in the rest of the 
calculations the in-plane component of the pump pseu- 
dospin is taken along the +x direction, i.e., 8 = 0). As 
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FIG. 4: Circular polarization degree of the driven mode versus 
the circular polarization degree of pump for different values 
of the longitudinal-transverse splitting Q for T=I00 K. Thus 
TE-TM splitting revives the bistability behavior even at high 
temperature values. At strong TE-TM splittings, i. e., Q, = 
1.5 meV, the circular polarizability diminishes due to strong 
mixing of a ± modes, and the hysteresis loop disappears. 



it was shown before in Fig. 2, in the absence of TE- 
TM splitting due to the dissipative nature of polariton- 
phonon interactions the hysteresis behavior is washed out 
at 100 K. On the other hand, if Q ^ 0, as in Fig. 4, 
the bistable behavior can be recovered and bistability 
phenomena can survive up to higher temperatures as 
compare to the case f2 = 0. However, instead of two- 
stepped hysteresis loop shown in Fig. [31 we observe only 
one-stepped behavior. This result suggests the transition 
from two-independent modes dynamics (the two modes 
are the mode with a + and the mode with o~ polariza- 
tions) to a single collective mode dynamics. In fact, at 
some critical value of the TE-TM splitting, (around 0.1 
meV in our parameter regime) , transition from high pop- 
ulation state to low population state of the one mode is 
always accompanied by the simultaneous transition from 
low population state to high population state of the other 
mode and crossover from the multistable behavior to the 
bistable occurs. 

If one increases the value of the TE-TM splitting 
field even further, polaritons would prefer to stay in 
quasi-lincarly polarized state due to strong mixing of cr* 
modes, even for pumping by the fully circularly polarized 
laser. This situation occurs at Q = 1.5 meV in Fig. 4, 
where polaritons become highly linearly polarized even 
at the values of p p = ±1 due to high value of the effec- 
tive magnetic field in the +x direction. The last term 
in Eq. (|2"5)) is responsible for this behavior. The parti- 
cles align their pseudospin parallel to the strong effective 
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magnetic field to minimize the total energy in the sys- 
tem. Meanwhile, the hysteresis behavior vanishes, and 
the difference between the backward and forward swap- 
pings disappears. 

IV. CONCLUSIONS 

In conclusion, we developed a formalism for the de- 
scription of the dissipative dynamics of an inhomoge- 
neous spinor polariton system in real space and time ac- 
counting for polariton-polariton interactions, polariton- 
phonon scattering and effect of the TE-TM effective mag- 
netic field. We applied our formalism to one-dimensional 
polariton condensate at different temperatures to inves- 
tigate the dynamics of the circular polarization of the 



system when it is driven by the external homogeneous 
laser pump. We showed that the polarization of the con- 
densate is highly sensitive not only to the history of the 
strength of the pump, but also to the phase of the ellip- 
tical polarization degree of this pump. In the presence of 
TE-TM field we observe the survival of this phenomena 
up to very high temperatures. 
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